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Abstract 

In search for a microscopic description of "droplet-like" properties for 
the Ising spin glass (single component order parameter, zero modes i.e. 
correlation functions vanishing at infinity) we reconsider the two-packet 
model of Bray and Moore, which is effectively Replica-Symmetric and 
enjoys zero modes but only up to one-loop. We show how an appropriate 
change in the limits of the basic parameters of the model (packet size 
and replica number) allows for a derivation of Ward-Takahashi (WT) 
identities, thus ensuring the existence of zero modes to all orders and 
opening the way for a Lagrangean formulation of a "droplet-like" field 
theory for the Ising spin glass. 

Spin glass theory has presented so far the schizophrenic aspect of two con- 
flicting approaches difficult to reconcile. The so-called mean-field like approach 
developed around Parisi P3|H| solution of the Sherrington-Kirkpatrick mean- 
field model uses standard field theory (mean-field, loops, renormalisation group, 
WT identities, ...) on a Lagrangean built with replicated fields. It is a "mi- 
croscopic" theory for the spin glass, in the same sense as the 4> 4 theory for 
the ferromagnet. The alternative droplet-like theory of Fisher and Huse [5] 
and Bray and Moore jS], despite abundant results, does not appeal to such a 
Lagrangean field theory starting point and is far removed from "microscopic" 
description. So there is a strong motivation to search better into the replicated 
4> 3 Lagrangean and examine whether the characteristic features for a "droplet" 
like theory could fit in. 

If one were to put up one, it would have the features of a theory with 
an effective Replica-Symmetric order parameter (a ferromagnetic in disguise). 
It would also possess zero modes to allow for an algebraic decrease at large 
distances of its correlations. But such features were indeed present, years before 
the birth of droplet theory, in an ansatz proposed by Bray and Moore 7 (BM), 
were replica symmetry was broken into "two-packets" , and restored at the very 
end. From their results, calculated at one-loop, it could be checked that both 
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features mentioned above (no RSB, zero modes) were present. However there 
was no guarantee that the zero modes would survive beyond one loop (indeed 
they would not and it was suggested that a generalization to multi-packet could 
help zero modes to remain massless). 

In this note we wish to take a new look at the "two-packets" RSB, as BM 
applied it in [7] to a spin glass Lagrangean with cubic coupling. We wish first 
to understand why this RSB did not give rise to WT identities that would have 
protected the zero modes to all orders. Indeed such identities have been derived 
8 in the framework of a Lagrangean field theory for systems with R steps of 
RSB (with in the end the Parisi limit B — » oo). One essential ingredient in the 
proof is the selection of "infinitesimal permutations" , in fact infinitesimal like 
1/R. In the two-packet theory, one packet has m replicas (a b c ...) and the 
other n — m (a f3 7 . . .). In the end n is set to zero (as it should for replicas) 
and m sent to infinity. Infinitesimal permutations are then easily selected (they 
are associated with transverse generators) and they go to zero with 1/m. So 
following the same steps as in jS], we identify below where the model fails to 
yield WT identities. Bringing the appropriate change allows then WT identities 
to be established, thereby giving life to a "droplet" Lagrangean field theory for 
the spin glass. 

1 Framework for WT identities: 

We start with a permutation invariant free-energy functional 

F{Qa,b} = F{Q pa ,pb} (1) 

A is a replica index that could have been denoted (i, a) i = 1 or 2, but is more 
economically replaced by a or a, the roman-greek notation of BM. The order 
parameter Qab is given by the stationarity condition on Q and as in , at 
mean-field level one has 
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On the other hand, from the definition of Legendre transform one has 1 

W {H AB } + F {Qab} = £ H ab Qab (3) 

(AB) 

where Hab is an external (unphysical) conjugate source, and hence 

dFjQcp} „ U] 
oQab 

1 For reasons of convenience we have changed the notation of BM with their Q3 — > Qi, 
Q2 — * Qoi Qi — > Q'2 
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yielding stationarity when the source is set to zero. 
The invariance under permutation then writes 

OF 

-rr^ — {Qcd + SQcd} = Hab + SH ab (5) 

OQAB 

where 

Qpa-pb = Qab + SQab (6) 

Now we have to make a choice for P (a choice amounting to use transverse 
generators). Just like in [5] we can divide packet one into equal bunches 
of p roman replicas and packet two into "~ T " equal bunches of greek replicas. 
The permutation chosen is for example, exchanging the first bunch of p roman 
replicas (ai b\ c\ ...) with the first bunch of p greek replicas (a\ (3\ 71 ...). In 
the following we keep a b c ... or a j3 7 ... for replicas that do no belong to the 
exchanged first bunches. 

Let us now look at the effect of such a permutation P by computing SQcd 
(or 5 Hab)- Clearly one has 

SQab = SQ a ji = 5Q aa = 5Q ai a t = (7a) 

5Q ai b t = -5Q ai f3 1 = Qs - Qi = SQo (7b) 

5Qa ia = SQaia = Ql ~ Qo = 5Q 2 (7c) 
SQ aai = SQ ai a = Ql ~ Qo = SQl (7d) 

and all these quantities are infinitesimal with 1/m. So that we are entitled 
to expand JSJ in SQcd and keep only the first term in its Taylor expansion, 
provided the resulting summation does not ruin the infinitesimality. Thus from 
101 we obtain, the formal WT identity 

E in ^Tn 6 Qcd = SH AB (8) 
OQabOQcd 

Note that this relationship has zero on its RHS for (AB) as in (|7a|l . or a non-zero 
RHS in SH for (AB) as in (|7EI ITcl lTcty . 

2 Some notations: 

To write out in detail eq. JSJ) we need to introduce some notation for 

8 2 F 

°^ m ab,cd (Q) 



OQabOQcd 
Noting the overlaps An B 



aHb=l (10a) 
a n (3 = 2 (10b) 
an« = (10c) 
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the matrix M AB ' CD will be identified by its overlaps AC\B and C C\D written 
as upper indices. To have a complete set of matrix elements we need also to 
specify how many maximal cross- overlaps we have: if AB ^ CD, 1 if A = C, 
or B = C, or A = L>, or B = D, and 2 if A = C, B = D or A = D, B = C. 
This closeness index we write as a heavy lower index. For example we have 

(11a) 

(lib) 

(11c) 

Alike for the 1 < — > 2 exchange in the upper indices. Also 



j^jab,ab 






Mi' 1 




M 1 ' 1 



M aa;aa = M 0;0 
M a»:bP = M 0;0 

The only ambiguity left is to distinguish M aa ' ab from M aa ' a P which we write 

(13a) 
(13b) 



M 



1(2) 



M aa;a0 = j^OjO 



exhibiting in parenthesis whether the repeated replica is roman (1) or greek (2). 

3 WT identity for AB = a\b : 

We carry out explicitly the J2cd summation. We have 

J Mai b; ai b + ^ M aib ' aiC + ^ M aib ' Mb + ^ M aib ' blC J (SQi) 

\ c bi bic / 

+ I ^M aib ' aia + J2M aib ' Ma J (SQ 2 ) 

\ a abi / 

+ ( ^M aifc;aifc + ^M aih;QlC j (SQi) 

\ ai aic / 

+ I M a ^ fi J {-SQ 2 ) = -8H X (14) 

Expliciting the summations and using the above notations (11)-(13) we get 
M 2 1;1 + (m - p - 1) Ml' 1 + (p - 1) Mj ;1 + (p - l) ( m - p _ l) M 1;1 1 5Q! 



(n - m - p) M^ } + (p-i)( n -m-p) M i;U <5Q 2 



1:0 



+ p(n-m-p) M 1 ' 2 = 8H 1 (15) 
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The first observation is that if we keep the limits n — > 0, followed by m — ► oo 
as taken in [7], the LHS of l|15|) contains factors going to infinity and it is no 
longer justified to replace eq. © by the first term in its Taylor expansion eq. 
|JH}. If we want to get rid of the terms in m in l|15f) we can choose 



(16) 



and with n — > first, and then /i — > oo. Note that a choice m = n Q /i, for 
a < 1/2 would sent Qo to infinity. As for the choice a > 1/2 it would imply 
Qi = Q2 = Qo = Q, leaving no room for the identities looked after. With the 
special choice a = 1/2, in the limit n = we have (2) replaced by 



h = Q2 

?o = Q4 



--Q 
Q_ 

2/i 2 



and (7) by 



SQo = 



SQi = 5Q 2 = - 7 ^ = SQ 



Q_ 

2fi 2 



(17a) 
(17b) 

(18a) 
(18b) 



Thus for /i large we have an infinitesimal transform. (Note that with (|18aH we 
did not bother to write terms in SQo that occur in eq. I|14|)l. 
With n = and with the choice of (|16|l we now get 



M 2 1;1 - 2M X 1;1 



M. 



1:1 



SQ + p 2 



-M. 



1:1 



2M 1;0 - Mo ;2 



6Q = 6H (19) 



Here p, the number of replicas in the exchanged bunch can be any finite integer, 
p > 1. So clearly, if we wish an unambiguous answer, it would be necessary that 
from other equations the coefficient of p 2 be set equal to zero. 



4 Related WT identities for AB = a x bi 

and AB = ab : 

Both identities have a vanishing RHS, aib\ leads to SHq (vanishing as in (|18al) ') 
and ab leads to zero (as in {7a|). The calculation follows exactly the same 
line as in the previous section. Spelling out the two equations obtained yields 
(i) the vanishing of the coefficient of p 2 in (|19l) , giving the diagonal component 
M Q ' in terms of the off-diagonal components Mq'°, A/q' 2 , (ii) the corresponding 
relationship for M x ' (see below). 



5 WT identities exhibited: 

With the above we can now express WT identities obtained for A n B = 1 : 
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Replicon for overlap 1 



r tt 

M 2 1;1 - 2M X 1;1 + M 1;1 = — (20a) 

and 

M 1;1 = 2M 1;0 - M 1;2 (20b) 
Mi ;1 = M^ + M o 1;0 - M 1;2 (20c) 

In exactly the same way one gets corresponding equations for A n B = 2 

M 2;2 - 2M 2;2 + M 2;2 = — (21a) 



M 2;2 = 2M 2;0 - M„ ;2 (21b) 

M 2;2 = M 2;0 + M 2;0 - M„ ;2 (21c) 

Finally taking iflB = 0, one gets 

M 2 °;° - (M^ + J^g)) + MT - % (22a) 



M ^ = i(M ^ + M ^ ) (22b) 
<S) = <° + i (m 2 * - M^°) (22c) 
M^ = M?° - i (m^ M?) (22d) 

6 Effect of a magnetic field: 

At mean-field level it is easily verified that the equations of motion that yield (2) 
are proportional to the equation giving the lowest eigenvalue (Replicon) of the 
Hessian. That is, the WT identity for the Replicon zero mode is trivially checked 
at the mean-field level. In the presence of an external magnetic field (distinct 
from the unphysical conjugate fields Hab) the equations of motion yielding 
the order parameters is unchanged but for an extra term H. This H cannot 
appear in the Hessian, a second derivative matrix, since in the Lagrangean it 
occurs in the linear term H Xmb 4>ab, with Qab — {<Pab)- Hence the presence 
of an external magnetic field suppresses the Goldstone modes and hence the 
transition, just like it occurs in O (N) systems. 
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7 A return on mean-field: 

Let us look back at eqs (2) which are only valid at the mean-field level. Actually, 
in our derivation, we only have used a milder form of eqs (18). To get the WT 
identities we only needed 

SQ a = (23a) 
SQ ~ 1/fj 2 (23b) 

for the limits 

n = (24a) 

■^<«1. (24b) 
A 1 

Is (23) valid to all orders beyond mean-field? This is easily checked to all orders 
in the paramagnetic phase. Going back to (|3I4I) . we have the order parameter 
Qab defined by 

Qab = dW {H} /dH AB (25) 

that is by 

QAB = Tr n D^CD <Pab exp I £{«£} + h cd<I>cd \ (26a) 

(CD) [ (CD) J 

M = [ II D(/)CD exp \ C ^> + J2 Hcd^cd \ (26b) 

(CD) [ (CD) J 

Here C {0} is the cubic BM Lagrangean, where the fields 4> A b (i) are coupled 
via |r J2i J2abc $ ab (*) fee (i) 4>ca («)• We have everywhere omitted the space 
(site) dependence since, in the end, the external source Hab (£) is always taken 
as Hab, site independent. 

Consider then the perturbation expansion of (25, 26) giving Qab as a power 
series of Hcd- If we choose Hcd = H, we then have Qab = Q and we write it 
as 

Q = f(H)H. (27) 

If we choose now H a b — H a p — H and H a p = Hq, one then has (no Hq 
dependence when n and m vanish) 

Qab = Qaf} = Q = f (H) H . (28) 

One also has (because W (H, Hq) can only be even in Hq) 

Q af 3 = Qo = g (H; H Q ) H Q (29) 
and from (|27|l when Ho = H, 

g(H;H) = f(H) . (30) 
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Hence for Hq — H ~ 1/^t 2 , we have 



/(ff) + «) (ff0 _ fl) 



and 



Q-Qo = 



f(H) + H 



dH 
dg(H;H) 



Ho 



dH 



(H-H )+O(H-H y 



That is, under the limits (24), one gets 

SQ ~ 6H 

thus justifying (23). 



(31) 
(32) 

(33) 



8 Comments and Conclusion: 

We have thus exhibited the Goldstone behavior for the three Replicon modes 
(20a, 21a, 22a), these modes acquiring a mass proportional to the (unphysical) 
conjugate field 8H, with in the end 8H — 0. A detailed examination of the 
Hessian matrix components shows |Hj besides that both the anomalous and 
longitudinal modes with zero overlap (as in (22)) also remain massless. 

Altogether we have 10 relationships (for 15 components). Note that with 
the five off-diagonal components one builds the seven diagonal terms Mq' 1 , Ml' 1 
with i = 1, 2, as in (20bc), (21bc), (22bcd). The last three diagonal terms 
M% (the one that contain the kinetic term for non-zero value of the momentum) 
complete the setup. 

To conclude we have given the right to exist to a spin glass theory whose 
starting point is formally identical to Bray and Moore two-packets theory but 
with the crucially different limits for the parameters n, m, namely with 



With the limit n = 0, fi — > oo we have then derived ten relationships between 
the fifteen components of the two-point (one particle irreducible) functions. 
Relationships between three-point functions could be derived in the same way 
all these relationships being non-perturbative. 

This new Lagrangean is a good candidate to describe "droplet" aspects of 
the Ising spin glass theory. It raises many questions and enjoys the following 
properties: 

(i) As in BM, the order parameter is, in the end, Replica Symmetric, a disguised 
fcrromagnet. 

(ii) Its associated free-energy is probably worse (lower) than the Parisi free- 
energy in the vicinity of the upper critical dimension. What would be the effect 
of dimension (that enters via loops) and would there be a critical dimension 
below which the "droplet" description would prevail is a crucial question to in- 
vestigate. 
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(iii) Its correlation functions enjoy several Goldstone modes. These modes be- 
come massive in the presence of an external magnetic field. They interact 
through cubic coupling. Thus their effective coupling cannot vanish in the infra- 
red like is the case for O (N) systems. It is thus expected that the 1/p 2 behavior 
of the Goldstone modes will develop anomalies. Large distance anomal behavior 
of correlation functions (computed for example in 6 — D dimension) will then 
have to be confronted with numerically obtained droplet exponents. 

The author is thanking E. Brezin, A. Crisanti and T. Temesvari for useful discussions. 
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2 With equation numbering having been messed up at the editing stage, the reader should 
rather consult the cond-mat version. 
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